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The inhomogeneous quark condensate, responsible for dynamical chiral symmetry breaking in cold 
nuclear matter, is studied by putting skyrmions onto the face-centered cubic crystal and treating the 
skyrmion matter as nuclear matter. By varying the crystal size, we explore the effect of density on 
the local structure of the quark-antiquark condensate. By endowing the light vector mesons p and uj 
with hidden local symmetry and incorporating a scalar meson as a dilaton of spontaneously broken 
scale symmetry, we uncover the intricate interplay of heavy mesons in the local structure of the 
quark condensate in dense baryonic matter described in terms of skyrmion crystal. It is found that 
the inhomogeneous quark density persists to as high a density as 4 times nuclear matter density. 
The difference between the result from the present approach and that from the chiral density wave 
ansatz is also discussed. 

PACS numbers: 11.30.Qc,11.30.Rd,12.39.Dc 

While the observation of the Higgs boson is considered 
to account for the masses of “elementary constituents” 
of visible matter, i.e., quarks and leptons, the bulk of the 
mass of the proton, the constituents of which are three 
light quarks, i.e., two up quarks and one down quark, 
remains more or less unexplained. The quark masses 
account for only < 1% of the proton mass. This is in 
marked contrast to the next scale hadron, the nucleus. 

The mass of a nucleus of mass number A is given al¬ 
most entirely, say, ~ 98%, by the sum of the masses of A 
nucleons. It is generally accepted, although not proven 
rigorously, that the nucleon mass arises from the non- 
perturbative dynamics of strong interactions encoded in 
QCD. What is involved is quark confinement and chiral 
symmetry spontaneously broken by the vacuum. It is 
not established but generally believed that the two are 
related. It is easier to address the latter in effective the¬ 
ories, so we will focus on it in what follows. 

In QCD, the order parameter of the chiral symmetry 
breaking is the pion decay constant /tt and, in terms of 
the fundamental QCD quantities, the quark condensates, 
typically the two-quark condensate (gg). The two-quark 
condensate in cold or hot dense matter has been stud¬ 
ied for a long time. In Refs. [il,i it was found that in 
neutron matter at the neutron star density, both the a 
(the isoscalar component of the chiral four-vector) and 
pion condensates exist and because of the chiral invari¬ 
ance of the system, these two condensates are linked to 
each other through a chiral rotation. Later, the approach 
put forward in Ref. [I| was investigated in more detail, 
and it was found that both condensates, a and the pion, 
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could depend on the position; that is, the condensates are 
inhomogeneous Q. In a phenomenological model with 
nucleons as explicit degrees of freedom which are put on 
a crystal lattice, Pandharipande and Smith found that in 
solid neutron matter, given an enhanced attractive tensor 
force, which is stronger than that given by standard one- 
pion exchange, the inhomogeneous pion condensate could 
take place y. Since then, the inhomogeneity of these 
condensates, now termed a chiral density wave(CDW), 
has been extensively studied in various different effective 
theory approaches modeling low-energy QCD; see, e.g., 
Refs. [1,0 and references therein. 

In this article, continuing our effort to understand 
hadronic matter at high densities, we approach this prob¬ 
lem by putting skyrmions, representing nucleons, on a 
face-centered cubic (FCC) crystal lattice. This has the 
advantage of enabling one to describe, on the same foot¬ 
ing, both the basic structure of the nucleon and the prop¬ 
erties of nuclear matter, thereby probing the modification 
of the vacuum in which the nucleons are propagating. 
By squeezing the crystal lattice, the density effect is sim¬ 
ulated and the properties of the condensates are modi¬ 
fied. The model we employ is the “generalized” skyrmion 
model that incorporates, in addition to the pion held car¬ 
rying topology, the vector mesons p and oj introduced as 
hidden local symmetry helds and the dilaton x associated 
with the broken scale symmetry of QCD. We call this 
the “dHLS skyrmion” to distinguish it from the Skyrme 
model with pion held only. It turns out that each of these 
helds plays an essential role in the structure of the nu¬ 
cleon as well as dense matter. In this article, we study the 
matter density effect and hadron resonance effect on the 
quark-antiquark condensates at low and high densities of 
nuclear matter. The details of what we obtained before 
are involved and given elsewhere (we refer the readers 
to Refs. 0 , [ 1 ]). We summarize the key points that are 
relevant to what we need for the discussions that follow. 
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1. There is a topology change from skyrmions to half- 
skyrmions at a density ni /2 with or without non- 
pion fields. It might involve no symmetry change 
of QCD and hence no order parameter like what 
happens in condensed matter physics The lo¬ 
cation of ni /2 depends, however, on the presence 
of the nonpion fields. Phenomenologically, it lies at 
ni /2 2no where uq is the normal nuclear matter 
density. 

2. The quark-antiquark condensate (qq) ^ 0 for n < 
ni /2 as is appropriate for the Nambu-Goldstone 
phase. At ni /2 the two-quark condensate, on av¬ 
erage, in a unit cell vanishes, so {qq) = 0 for 
n>nii 2 (here and in what follows, () stands for the 
space average in a unit crystal cell). However, the 
pion decay constant does not vanish, with massive 
hadrons excited, until a higher density ric > nij 2 i 
with Tic being the density at which the chiral sym¬ 
metry is restored. What happens is that the two- 
quark condensate is not locally zero but vanishes, 
on average, in the unit cell. Thus chiral symmetry 
is not actually restored. This means that the two- 
quark condensate is not an order parameter for chi¬ 
ral symmetry in this model. There must be differ¬ 
ent order parameters for chiral symmetry, perhaps 
of multiquark nature, with the Gell-Mann-Oakes- 
Renner relation still holding in the phase with pos¬ 
sible mutiquark condensates as an order parameter. 
The locally nonvanishing two-quark condensate im¬ 
mediately suggests the existence of an inhomogene¬ 
ity in the half-skyrmion phase. 

3. The nucleon mass can be characterized by two com¬ 
ponents asTTiN = mo + M{{qq)), where M{{qq)) —>■ 
0 as {qq) —0 and mo 0. Here, mo resembles the 
chiral-invariant mass term in the parity-doubled 
nucleon model Furthermore, as will be men¬ 
tioned later, parity doubling does take place in the 
hadron spectrum in a cold or hot dense medium. 
However, in the dHLS model, the mo term which is 
related to the chiral doubling structure of nucleons 
does not figure explicitly in the Lagrangian, so it 
most likely reflects an emerging symmetry induced 
by the medium. 

4. There is an intricate interplay between the isovec¬ 
tor vector meson p, the isoscalar vector meson w, 
and the scalar (dilaton) meson in dense baryonic 
matter. Suppose one generalizes the model used 
here by putting the infinite tower of isovector vector 
mesons p, p' etc. and oi, a( etc. in the skyrmion 
model while dropping all isoscalar vector mesons 
and scalar mesons. Then one can arrive in flat 
space at many-nucleon systems - and infinite mat¬ 
ter - that are described by BPS skyrmions with 
vanishing binding energies M- This is close to 
what is observed in experiments for binding ener¬ 
gies of medium and heavy nuclei. What one notes 


is that each member of the tower of the isovector 
mesons tends to move the system closer to matter 
of BPS skyrmions. In the model used in this arti¬ 
cle, however, this tendency with the p is spoiled by 
the isoscalar vector meson wQ 

In an effective model of QCD with linear chiral sym¬ 
metry, the chiral four-vectors in the meson fields M are 
decomposed as 

M{x) = a{x) + iTana{x), with a =1,2,3, (1) 

with r“ as Pauli matrices. In the Nambu-Goldstone 
phase, the chiral four-vectors are constrained by (cr^ + 
7>'a)QCD = /I (here and in what follows, ()qcd stands 
for the medium-free QCD vacuum) in the QCD vacuum. 
Moreover, due to the parity invariance of the QCD vac¬ 
uum, one normally has ((t)qcd = U oc (gg)QCD and 
(7ra)QCD OC (glTa75g)QCD = 0. However, in a medium, 
due to the interaction between the mesons and baryonic 
matter, both (a) and (tt^) could be nonzero and position 
dependent [H-Q. 

In the nonlinear realization of chiral symmetry, the 
two-quark condensate is normalized as {U)qcd = 
{M)qcd/Jtt = 1 in the QCD vacuum with U = 
exp(i7r//^). However, in the nuclear matter medium, 
there can be deviation due to the interaction between 
the pion and the matter medium, and the value of the 
deviation from I accounts for the deviation of the two- 
quark condensate from its vacuum value. 

To explore the density dependence of the quark- 
antiquark condensate in a theory with the nonlinear real¬ 
ization of chiral symmetry, it is convenient to decompose 
the field U{x) as in JT]), 

U{x) = (l)oix) + iTaipaix), with 0=1,2, 3. (2) 

In the ground state of the system, 4io{x) accounts for the 
isoscalar quark condensate while (paix) accounts for the 
isovector quark condensate which, in the QCD vacuum, 
have values 4>o{x) = 1 and (j)a{x) = 0, but in the nu¬ 
clear matter medium, they might have 4>q{x) ^ 1 and 
^a{x) 7 ^ 0 due to the matter effect. And, because of the 
cluster structure of nuclear matter, the quark conden¬ 
sates are expected to be functions of space coordinates 
X, i.e., inhomogeneous. 

In the present work, to simulate the nuclear matter 
environment, among all the approaches to nuclear mat¬ 
ter, we adopt the skyrmion matter approach and regard 


^ The model with the infinite tower of isovector vector mesons 
arises in the large Nc and large’t Hooft constant (A) limit by di¬ 
mensional deconstruction from 5D h olog raphic dual QCD model 
that arises from string theory [TH . Il2l . With next-order 1/A 
corrections, the isoscalar vector mesons coming from the Chern- 
Simons term and the warping of the metric, together, destroy 
the BPS structure, so one loses the nice agreement with nature. 
How this pro blem can be resolved is another matter discussed 
elsewhere fU^ . 
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the skyrmion matter as nuclear matter [l^ . [l5l | . Since in 
each crystal cell, the baryon number is fixed, the density 
effect is simulated by charigingthe crystal size. Here, 
we adopt the FCC crystal [la, [T^ for which the nuclear 
matter density n and crystal size L are related through 
the relation n = 4/(2L)^. As stated at the beginning, 
in the skyrmion crystal approach, one can treat both the 
nuclear matter and medium modified hadron properties 
in a unified way [l8l |. 

To simulate nuclear matter from FCC crystal, it is con¬ 
venient to define the unnormalized quantities pin, 
which are related to the quantities in the decomposition 
([2]) through 

=, (a,^ = 0,1,2,3). (3) 

\/EJ.„(0s) 

In the crystal lattice, the unnormalized fields have the 
Fourier series expansions 

Mx,y,z) = E Pabc cos(a7rx/T) cos(67ry/T) cos(c7rz/L), 

a,6,c 

E ahki siii{h7Tx /L) cos{k'Ky /L) cos{1’kz/L)^ 

h,k,l 

^2{x,y,z) = E ahki cos{l7Tx/L) sin(ft,7ry/L) cos{kTrz/L), 

h,k,l 

h{x,y,z) = E ahki cos{k'Kx/L) cos{l7Ty/L) sin^hirz/L). 
h,kd 


By varying the Fourier coefficients a and /3, one obtains 
the minimal energy, interpreted as the ground state of 
nuclear matter, for a specified crystal size L. The inho¬ 
mogeneous quark condensates with respect to their val¬ 
ues in vacuum can then be calculated by substituting the 
values of a and /3 thus determined into Eq. ([3]). In addi¬ 
tion, by varying the crystal size L, the density effect on 
the quark condensate can be accessed. Note that due to 
the FCC structure and the arrangement of the nearest- 
neighbor skyrmions to yield the strongest attractive in¬ 
teraction, the modes appearing in the above equation are 
not independent [T^. 

To explore the tt, p, and w meson effect on the inho¬ 
mogeneous quark condensate, we use the skyrmion model 
obtained from the hidden local symmetry approach given 
up to the next-to-leading order in the chiral expansion 
including the homogeneous Wess-Zumino term. The de¬ 
grees of freedom involved are then the pion, the lowest- 
lying p and w as employed in to explore the struc¬ 

ture of the single baryon as well as baryonic matter. Ex¬ 
plicitly, we consider the models HLS(7r), HLS(7r,/A and 
HLS(7r,p,w) which are explicitly defined in Refs. [3,Hi' 
In addition to the hidden local symmetry degrees of free¬ 
dom, to explore the light scalar meson (i.e., dilaton) effect 
on nuclear matten we apply the model dHLS-II(7r, p, lu) 


defined in Ref. |8| in which the chiral restoration can 
be realized at some high density ric- The effective La- 
grangian of model dHLS-II(7r, p, w) is 


r — /-dHLSII 

E-dHLSII — -^(2) 




HLS 

(4) 


’HLS 


-'dilaton 5 


( 4 ) 


where and are the O(p^) and homogeneous 

Wess-Zumino terms of the HLS Lagrangian, respectively, 
and 

.^dilaton = 

771 ^ 

_ X^X 

4 

+ o/J (7) Tr[d||„6j|su(2) + 

+ -^kin, (5) 

with the subindex SU(2) being the SU(2) subgroup of 
the hidden local symmetry U(2) and £kin being the ki¬ 
netic term of the dilaton field. In the Lagrangian (|4|), 
when the dilaton field is replaced by its vacuum expecta¬ 
tion value, (x) = f^, the model dHLS-H(7r, p, w) reduces 
to HLS(7r,p,w). After switching off the homogeneous 
Wess-Zumino term, one arrives at HLS(7r,p) and, after 
integrating out the p meson from HLS(7r, p) we obtain the 
model HLS(7r). Eor other details of the relevant models, 
we refer to Refs. 0, 

From Eq. m we see that, in the skyrmion crystal ap¬ 
proach, all the quark condensates are functions of the 
three-dimensional coordinates. In Fig. [T] we illustrate the 
crystal size dependence of the quark-antiquark conden¬ 
sates (j)o{x) oc {qq(x))* and (j)i(x) oc {qT^j 5 q{x))* at the 
z = 0 plane using the model HLS(7r,p, w) in which all 
the TT, p, and lu meson contributions are included and a 
self-consistent power-counting mechanism exists. These 
plots tell us that in both the skyrmion and half-skyrmion 
phases, the chiral symmetry is locally broken since the 
quark condensates could be locally nonzero. 

Let us first look at the first row of Fig. [TJ From 
Fig.dja) one sees that the large area is occupied by white 
space so that at this density the space average of (/)o is 
approximately its value in the QCD vacuum, i.e., ^ 1, 
which agrees with the result of Ref. Q. By comparing 
Figs. Hla) and [TJb) , we see that a larger space is occu¬ 
pied by negative values of (po, indicating that the space 
average of po becomes smaller as density increases. The 
comparison of Fig. [TJb) - which is in the skyrmion phase 
- and Fig. [TJc) - which is in the half-skyrmion phase - 
indicates that the changeover from the skyrmion phase 
to the half-skyrmion phase could be smooth. Finally, the 
close resemblance of Figs. [He) and [TJd) reflects the fact 
that even in the half-skyrmion phase the chiral symme¬ 
try is still locally broken, though the condensate is zero 
averaged in the unit cell and the order parameter for chi¬ 
ral symmetry breaking in medium /* is nearly constant, 
independent of density. 

Next we consider the second row of Fig. [T] which 
describes the density dependence of the inhomoge¬ 
neous parity-odd isospin nonsinglet quark condensate 
(gr^ysg)*. We first find that this inhomogeneous conden- 
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(a) L = 2.5 fm 



(b) L = 1.5 fm 


(c) L = 1.4 fm 
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FIG. 1. (Color online) The density effect on (j>o{x,y,0) (first row) and (j>i{x,y,0) (second row). The half-skyrmion phase 
appears at L = 1.45 fm. 


sate is an odd function of the coordinate x because in the 
skyrmion crystal approach the flipping of the coordinate 
Xi is accompanied by the flipping of the corresponding 
component of (j)i and, because of this property, we con¬ 
clude that the space average of the quantity (jji is zero. 
The figures in this row also show the density dependence 
of the component (j)i in the skyrmion phase and its nearly 
density independence in the half-skyrmion phase. 

We plot in Fig. [2] the effect of the (jj meson and the 
scalar (dilaton) meson on the quark condensates in the 
x-y plane. In this plot, we take the typical crystal size 
L = 1.5 fm at which the system is in the skyrmion phase 
for all four models. Compared to the model with a pion 
only, HLS(7r), we find that the plot from HLS(7r,p) is 
occupied by less blue (soliton matter) space, which indi¬ 
cates that the quark condensates appear in a smaller area 
due to the attractive force from the p meson. Similarly, 
the subtle difference between Figs. IKc) and [2Id) tells 
us that the inclusion of the dilaton field associated with 
the dynamical breaking of the scale symmetry slightly 
shrinks the area of the inhomogeneous quark condensates 
due to the attractive force from the dilaton. However, 
the dramatic difference between the results in Figs. [2](b) 
and me) result from HLS(7r, p) and HLS(7r, p, w) shows 
the crucial effect on the quark condensates from the lj 
meson. The repulsive force coming from the lj exchange 
causes the inhomogeneous quark condensates to expand. 
This may explain the deviation from the BPS structure 
caused by the presence of the uj meson noted in item S) 

As mentioned, in our model, the chiral condensate is 
locally broken in both the skyrmion phase, and half- 
skyrmion phase although in the latter phase it is re¬ 
stored globally. There is a basic difference between this 
skyrmion prediction and the CDW phase discussed in 
Ref. [^. First of all, in the quark-antiquark conden¬ 
sate is the bona-fide order parameter indicating the onset 


of the Wigner phase, whereas it is not in the skyrmion 
description, where it is more like a “pseudo-gap”. Sec¬ 
ond, the “dilaton” field x in @ is what corresponds to 
the isoscalar tetraquark state or a pion-pion resonance, 
whereas our dilaton is /o(500) which can be associated 
with a pseudo-Goldstone boson in the vicinity of a QCD 
IR fixed point [2lj |. In our scheme, there is no low-lying 
two-quark configuration fluctuation at low density cor¬ 
responding to the fourth component of the chiral four- 
vector in the linear realization of chiral symmetry. 

It is, however, interesting to make a comparison be¬ 
tween the two models. To do so, we take the Dautry- 
Nyman ansatz Q as done in 

(j)o{x) = cos(2/x); (/)3 = sin{2fx). (6) 

In this form, the parameter / is an indicator for the inho¬ 
mogeneous quark condensates. It was found in @ that, 
at low chemical potential, / = 0; i.e., there are no inho¬ 
mogeneous quark condensates. However, at high chemi¬ 
cal potential, / becomes a nearly nonzero constant which 
indicates the onset of the inhomogeneous quark conden¬ 
sates. In terms of the chiral field U, the indicator / in 
Q is simply 

f =-'-Tr (r^U^^d^Uo) . (7) 

However, in the skyrmion crystal approach, the situa¬ 
tion becomes somewhat complicated because of both the 
three-space dependence and the three-isospin dependence 
locked into the hedgehog. Here the relevant quantity is 

f^j{x,y,z) = -^Tr (riU^djUo^ 

~ ^ “b ^iab4^a^j4^b\ ; (S) 

which is a 3 X 3 matrix in the spin-isospin space, and each 
matrix element depends on the coordinates x,y, and z. 
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(a)HLS(7r) 


(b)HLS(7r, p) 


(c)HLS(7r, p,ui) 


(d)dHLS-II(7r,p,a;,x) 


FIG. 2. (Color online) The p, u) and x meson effects on (j>o{x,y,0) (first row) and <j)i{x,y,0) (second row) in the skyrmion 
phase at L = 1.5 fm. 


An explicit numerical simulation shows that in the FCC 
crystal, the quantities fij{x,y,z) vanish in neither the 
skyrmion phase nor the half-skyrmion phase, so the in¬ 
homogeneous quark condensate exists locally in both the 
skyrmion phase and the half-skyrmion phase. The inho¬ 
mogeneity with / 7 ^ 0 in the skyrmion phase could be 
due to the fact that at low density the classical skyrmion 
configuration in crystalline form is not in the Fermi liq¬ 
uid state. The homogeneity may be recovered by collec¬ 
tive quantization to restore spin and isospin. The half- 
skyrmion phase at n > ni/ 2 , however, is likely non-Fermi 
liquid in a crystalline form. 

In the CDW calculation studied in Ref. the indi¬ 
cator of the homogeneous quark condensate is the quan¬ 
tity 4> defined there related to the pion decay constant 
via /,r = with Z being the wave-function renor¬ 

malization constant of the pseudoscalar mesons arising 
from removing the tt-ui mixing terms. It was found that 
in the low chemical potential region, the homogeneous 
quark condensate (j) decreases slightly with the increase 
of the chemical potential, while in the high chemical po¬ 
tential region, it remains as a constant with about 25% 
of the vacuum. However, in the FCC skyrmion crys¬ 
tal approach, in all the models considered in the present 
work, the homogeneous quark condensate accounted for 
by the space-averaged a decreases with increasing density 


in the skyrmion phase and vanishes in the half-skyrmion 
phase 0 , 11 ]. 

Although the inhomogeneity in the quark-antiquark 
condensate persists in both the skyrmion phase and the 
half-skymrion phase with no obvious phase transitions, 
there is a drastic change in nuclear physics across the 
density ni /2 [111. As mentioned, parity doubling seems 
to emerge in the hadronic spectrum suggested in the 
parity-doubliM effective Lagrangian model and lat¬ 
tice studies [24| . Indeed, it was shown very recently in 
Ref. [Ill by coupling the heavy-light mesons with chiral- 
partner structure to the leading-derivative order HLS(7r), 
that the masses of the chiral partners of the heavy-light 
mesons, which are split in the skyrmion matter, become 
degenerate in the half-skymion phase. It would be in¬ 
teresting to study similar phenomena in the light meson 
sector. 
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